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ABSTRACT 

The last few months have seen the discovery of a second Martian Trojan 
(1998 VF31), as well as two further possible candidates (1998 QH56 and 1998 
SD4). Together with the previously discovered Martian satellite 5261 Eureka, 
these are the only known possible solar system Trojan asteroids not associated 
with Jupiter. Here, maps of the locations of the stable Trojan trajectories 
of Mars are presented. These are constructed by integrating an ensemble of 
in-plane and inclined orbits in the vicinity of the Martian Lagrange points 
for between 25 million and 60 million years. The survivors occupy a band of 
inclinations between 15° and 40° and longitudes between 240° and 330° at the 
L5 Lagrange point. Around the L4 point, stable Trojans inhabit two bands 
of inchnations (15° < z < 30° and 32° < i < 40°) with longitudes restricted 
between 25° and 120°. Both 5261 Eureka and 1998 VF31 lie deep within one 
of the stable zones, which suggests they may be of primordial origin. Around 
Mars, the number of such undiscovered primordial objects with sizes greater 
than 1 km may be as high as ~ 50. The two candidates 1998 QH56 and 1998 
SD4 are not presently on Trojan orbits and will enter the sphere of influence of 
Mars within half a million years. 

Subject headings: Solar System: general - planets and satellites: general - minor 
planets, asteroids 
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1. INTRODUCTION 

The Lagrange points L4 and L5 are stable in the restricted three body problem (e.g., 
Danby 1988). However, the long-term survival of Trojans around the Lagrange points of 
the planets in the presence of perturbations from the remainder of the Solar System is a 
difficult and still unsolved problem (e.g., Erdi 1997). Jovian Trojan asteroids have been 
known since the early years of this century, while a number of Saturnian moons (e.g., Dione 
and Helene, Tethys and Calypso, Tethys and Telesto) also form Trojan configurations with 
their parent planet. However, whether there exist Trojan-like bodies associated with the 
other planets has been the matter of both observational activity (e.g., Tombaugh 1961, 
Kowal 1971) and intense theoretical speculation (e.g., Weissman & Wetherill 1974, Mikkola 
& Innanen 1990, 1992). The answer to this problem came in 1990, with the discovery of 
5261 Eureka, the first Trojan around Mars (see Mikkola et al. 1994 for details). The last 
few months of 1998 have seen further remarkable progress with the discovery of one certain 
Martian Trojan, namely 1998 VF31, as well as two further candidates, namely 1998 SD4 
and 1998 QH56 (see the Minor Planet Electronic Circulars 1998-W04, 1998-R02, 1998-S20 
and the Minor Planet Circular 33085). The suggestion that 1998 QH56 and 1998 VF31 
might be Martian Trojans was first made by G.V. Williams. 

These recent discoveries raise very directly the following questions. Are there any more 
Martian Trojans? If so, where should the observational effort be concentrated? Of course, 
the first question can only be answered at the telescope, but the resolution of the second 
question is provided in this Letter. By integrating numerically an ensemble of inclined 
and in-plane orbits in the vicinity of the Martian Lagrange points for 25 and 60 million 
years respectively, the stable regimes are mapped out. On re-simulating and sampling 
the ensemble of stable orbits, the probability density of Martian Trojans as a function of 
longitude and inclination can be readily obtained. If a comparatively puny body such as 
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Mars possesses Trojans, the existence of such objects around the larger terrestrial planets 
also merits very serious attention. There are Trojan orbits associated with Venus and the 
Earth that survive for tens of millions of years (e.g., Tabachnik & Evans 1998). If objects 
populating such orbits exist, they must be small else they would have been found by now. 



2. MARTIAN TROJANS 

Saha & Tremaine (1992, 1994) have taken the symplectic integrators developed by 
Wisdom & Holman (1991) and added individual planetary timesteps to provide a fast code 
that it is tailor-made for long numerical integrations of low eccentricity orbits in a nearly 
Keplerian force field. In our simulations, the model of the Solar System consists of the eight 
planets from Mercury to Neptune, together with test particles starting near the Lagrange 
points. The effect of Pluto on the evolution of Martian Trojans is quite negligible. Of course, 
the Trojan test particles are perturbed by the Sun and planets but do not themselves exert 
any gravitational forces. The initial positions and velocities of the planets, as well as their 
masses, are provided by the JPL Planetary and Lunar Ephemerides DE405 and the starting 
epoch is JD 2440400.5 (28 June 1969). All our simulations include the most important 
post-Newtonian corrections, as well as the effects of the Moon. Individual timesteps are 
invaluable for this work, as orbital periods are much smaller in the Inner Solar System than 
the Outer. For all the computations described in this Letter, the timestep for Mercury 
is 14.27 days. The timesteps of the planets are in the ratio 1 : 2 : 2 : 4 : 8 : 8 : 64 : 64 
for Mercury moving outwards, so that Neptune has a timestep of 2.5 years. The Trojan 
particles all have the same timestep as Mercury. These values were chosen after some 
experimentation to ensure the relative energy error has a peak amplitude of ~ 10~^ over the 
tens of million year integration timespans. After each timestep, the Trojan test particles 
are examined to see whether their orbits have become hyperbolic or if they have entered 
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the planet's sphere of influence (defined as Tg = OpMp/^ where Op and Mp are the semimajor 
axis and mass of the planet). If so, they are terminated. In methodology, our calculations 
are very similar to the magisterial work on the Trojan problem for the four giant planets by 
Holman & Wisdom (1993). The earher calculations of Mikkola & Innanen (1994, 1995) on 
the Trojans of Mars for timespans of between tens of thousands and 6 million years have 
also proved influential. Our integrations of Trojan orbits are pursued for durations ranging 
from 25 to 60 million years, the longest integration periods currently available. Nonetheless, 
the orbits have been followed for only a tiny fraction of the age of the Solar System (~ 4.5 
Gigayears), so it is wise to remain a httle cautious about our results. 

Figure 1 shows the results of our first experiment. Here, the orbits of 1080 Trojan test 
particles around Mars are integrated for 25 million years. The initial inclinations of the test 
particles (with respect to the plane of Mars' orbit) are spaced every 2° from 0° to 90° and 
the initial longitudes (again with respect to Mars) are spaced every 15° from 0° to 360°. 
The starting semimajor axes and the eccentricities of the Trojans are the same as the parent 
planet. Only the test particles surviving till the end of the 25 million year integration are 
marked on the Figure. The survivors occupy a band of inclinations between 10° and 40° and 
longitudes between 30° and 120° (the L4 Lagrange point) or 240° and 330° (the L5 point). 
On the basis of 4 million year timespan integrations, Mikkola & Innanen (1994) claim that 
stable Martian Trojans have inclinations between 15° and 30° and between 32° and 44° with 
respect to Jupiter's orbit. Our longer integrations seem to suggest a more complex picture. 
Mikkola & Innanen's instabihty strip between 30° and 32° can be detected in Figure 1, but 
only for objects near L4 with initial longitudes < 60°. In particular, this instability strip 
does not exist around L5 and here Trojans with starting inclinations 30° < i < 32° seem to 
be stable - as is also evidenced by the recent discovery of 1998 VF31. Marked on the figure 
are the instantaneous positions of the two certain Martian Trojans, namely 5261 Eureka 
(marked as a red circle) and 1998 VF31 (a green circle), as well as the two candidates 1998 
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QH56 (a blue circle) and 1998 SD4 (a yellow circle). It is delightful to see that the two 
securely established Trojans lie within the stable zone, which was computed by Tabachnik 
& Evans (1998) before the discovery of 1998 VF31. In fact, they live deep within the heart 
of the zone, suggesting that they may even be primordial. The two candidates (1998 QH56 
and 1998 SD4) lie closer to the rim. Let us finally note that Trojans starting off in or near 
the plane of Mars' orbit are unstable. This has been confirmed by an extensive survey of 
in-plane Martian Trojans. On integrating 792 test particles with vanishing inclination but 
with a range of longitudes and semimajor axes, we found that all are unstable on timescales 
of 60 miUion years. Martian Trojans with low inclinations are not expected. 

It is useful to an observer hoping to discover further Trojans to provide plots of the 
probability density. Accordingly, let us re-simulate the stable zones with much greater 
resolution. This is accomplished by placing a total of 746 test particles every 1° in initial 
inclination and every 5° in initial longitude so as to span completely the stable regions. 
This ensemble of orbits is then integrated and the orbital elements are sampled every 2.5 
years to provide the plots displayed in Figure 2. The upper panel shows the meshed surface 
of the probability density as a function of both inclination to the invariable plane and 
longitude with respect to the planet. The asymmetry between the two Lagrange points is 
evident. The lower panels show the projections of the meshed surface onto the principal 
planes - in particular, for the inclination plot, we have shown the contribution at each 
Lagrange point separately. There are a number of interesting conclusions to be drawn from 
the plots. First, as shown by the dotted fine, the probability density is bimodal at L4. It 
possesses a flattish maximum at inclinations between 15° and 30° and then falls sharply, 
before rising to a second maximum at 36°. At L5, all inclinations between 15° and 40° carry 
a significant probability, though the smaller inclinations in this band are most favored. It is 
within these inclination windows that the observational effort should be most concentrated. 
Second, the probability density is peaked at longitudes of ~ 60° (L4) and ~ 300° (I/5). The 
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most likely place to observe one of these Trojans is indeed at the classical locations of the 
Lagrange points. This is not intuitively obvious, as any individual Trojan is most likely to 
be seen at the turning points of its longitudinal libration. There are two reasons why this 
effect is not evident in our probability density plots. First, our figures refer to an ensemble 
of Trojans uniformly populating the stable zone. So, the shape of the stable zone also plays 
an important role in controlling the position of the maximum of the probability density. 
Second, the positions of the Lagrange points themselves are oscillating and so the turning 
points of the longitudinal libration do not occur at the same locations, thus smearing out 
the enhancement effect. 

Table 1 lists the orbital elements of the two secure Martian Trojans and the two 
candidates, as recorded by the Minor Planet Center. Prom the instantaneous elements, it 
is straightforward to simulate the trajectories of the objects. Pigure 3 shows the orbits 
plotted in the plane of longitude (with respect to Mars) versus semimajor axis. As the 
figures illustrate, both 5261 Eureka and 1998 VP31 are stable and maintain their tadpole 
character (see e.g., Garfinkel 1977) for durations of 50 million years. Based on preliminary 
orbital elements, Mikkola et al. (1994) integrated the orbit of 5261 Eureka and found 
that its longitudinal libration was large, quoting 297° ± 26° as the typical range in the 
longitudinal angle. Our orbit of 5261 Eureka, based on the latest orbital elements, seems to 
show a smaller libration of 285 — 314°. The remaining two objects that have been suggested 
as Martian Trojans, 1998 QH56 and 1998 SD4, both enter the sphere of infiuence of Mars - 
in the former case after ~ 500 000 years, in the latter case after ~ 100 000 years. Although 
the orbits are Mars crossing, their eccentricities remain low and their inclinations oscillate 
tightly about mean values until the Mars' sphere of influence is entered. It is possible that 
these objects were once Trojans and have been ejected from the stable zones, a possibility 
that receives some support from their locations in Pigure 1 at the fringes of the stable 
zones. Of course, another possibility is that they are ejected asteroids from the Main Belt. 
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The fact that both confirmed Martian Trojans he deep within the stable zones 
in Figure 1 suggests that these objects may be primordial. If so, we can get a crude 
estimate of possible numbers by extrapolation from the number of Main Belt asteroids 
(c.f., Holman 1997, Evans & Tabachnik 1999). The number of Main Belt asteroids A^mb is 
A^MB ^ Smb^mb/, where Amb is the area of the Main Belt, Emb is the surface density of 
the proto-planetary disk and / is the fraction of primordial objects that survive ejection 
(which we assume to be a universal constant). Let us take the Main Belt to be centered on 
2.75 AU with a width of 1.5 AU. The belt of Martian Trojans is centered on 1.52 AU and 
has a width of < 0.0025 AU. If the primordial surface density falls off inversely proportional 
to distance, then the number of Martian Trojans Nut is 



The number of known Main Belt asteroids with diameters > 1 km is > 40000, which 
suggests that the number of Martian Trojans is ^ 50. 



Motivated by the recent discovery of a new Mars Trojans (1998 VF31) as well as further 
possible candidates (1998 QH56, 1998 SD4), this paper has provided maps of the stable 
zones for Martian Trojans and estimates of the numbers of undiscovered objects. For Mars, 
the observational effort should be concentrated at inchnations satisfying 15° < i < 30° 
and 32° < i < 40° for the L4 Lagrange point and between 15° and 40° for L5. These 
are the spots where the probability density is significant (see Figure 2), though the lower 
inclinations in these bands are slightly more favored than the higher. Trojans in or close 
the orbital plane of Mars are unstable. Crude estimates suggest there may be as many as 
~ 50 undiscovered Martian Trojans with sizes ^ 1 km . The orbits of 5261 Eureka and 1998 
VF31 remain Trojan-like for durations of at least 50 million years. The other candidates. 
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CONCLUSIONS 
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1998 QH56 and 1998 SD4, are not currently Trojans, though it is conceivable that they 
may once have been. Both objects will probably enter the sphere of influence of Mars after 
^0.5 million years. 
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Fig. 1. — This figure shows the stabihty zones of the inchned Trojans of Mars. The horizontal 
axis marks the longitude measured from Mars and the vertical axis the inclination with 
respect to Mars of the starting positions of test particles. At outset, the array of particles 
has inclinations spaced every 2° and longitudes spaced every 15°. The initial semimajor axes 
and eccentricities of the Trojans are the same as Mars. Only the particles surviving till the 
end of the 25 million year integration are marked on the figure, which provides a map of the 
stable regions. All the objects starting in- plane do not persist and only the inclined Trojans 
are stable. Also marked on the figure are the instantaneous positions of the two Martian 
Trojans, namely 5261 Eureka (marked as a red circle) and 1998 VF31 (a green circle), as 
well as the asteroids 1998 QH56 (a blue circle) and 1998 SD4 (a yellow circle). 

Fig. 2. — These figures show the most likely places to observe new Martian Trojans. They 
display the two-dimensional probability density as a function of the inclination with respect 
to the invariable plane and the longitude with respect to Mars (upper panel) together with 
the projections onto the principal planes (lower panels). The figures are constructed by 
re-simulating the stable regions displayed in Figure 1 at much greater resolution. 746 
test particles are placed every 1° in inclination and every 5° in longitude so as to span 
the stable region and the trajectories are sampled every 2.5 years for 50 000 years. The 
overall normalisation of the probability density is arbitrary. In the inclination plots, the 
contributions from the Lagrange points are separated - broken lines refer to L4 and unbroken 
lines refer to L5. 

Fig. 3. — Plots of the longitude versus semimajor axis are shown for the orbits of 5261 Eureka 
(upper panel) 1998 VF31 (lower panel). The orbits are integrated for 50 million years and 
are sampled every 10 000 years. 
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Asteroid 


a 


e 


i [deg] 


A [deg] 


H 


Diameter 


5261 Eureka 


1.5235 


0.0647 


20.2802 


301.4 


16.1 


2-4 km 


1998 VF31 


1.5241 


0.1005 


31.2966 


290.9 


17.1 


1-2 km 


1998 QH56 


1.5506 


0.0307 


32.2222 


258.5 


17.9 


1-1.5 km 


1998 SD4 


1.5149 


0.1254 


13.6725 


245.5 


18.7 


0.5-1.2 km 



Table 1: This table lists some of the properties of the two definite Martian Trojans, as well 
as two suggested candidates. These include the instantaneous semimajor axis a, eccentricity 
e, inclination from the J2000 plane i and longitude measured from Mars A. The epoch is JD 
2451200.5 (22 January 1999). The magnitude H and the approximate diameter of the object 
(inferred using albedos of 0.05 — 0.25) are also given. Most of this information is abstracted 
from Minor Planet Circulars 30250 and 33085 (Eureka and 1998 QH56) and Minor Planet 
Electronic Circular 1998-W04 and 1998-S20 (1998 VF31 and 1998 SD4). 
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Asteroid 


Ao 


Ae 


Ai [deg] 


AA [deg] 


Lag. point 


D [deg] 


5261 Eureka 


5.955 X 10-^ 


0.08117 


4.8874 


29.47 


7- 14.3 
-^5-15.2 


12.0862 


1998 VF31 


1.386 X 10-3 


0.1128 


4.9722 


69.38 


T +30.1 
-'^5-39.2 


47.523 



Table 2: This table lists some of the properties of the orbits of the two confirmed Martian 
Trojans, inferred from numerical integrations. The table gives the maximum variation during 
the entirety of the 50 million year integration timespan in the semimajor axis Aa, in the 
eccentricity Ae, in the inclination Ai and in the longitude measured from Mars AA. Both 
Trojans oscillate around the Lagrange point L5 and the superscript and subscript indicate 
the extent of the angular libration. Part of this includes the oscillation of the Lagrange 
point itself, so the final column D is the peak to peak angular libration measured from the 
Lagrange point. 
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